Comparison of methods for analysis of selective genotyping survival data

Survival traits and selective genotyping datasets are typically not normally distributed, thus common models used to identify QTL may not be statistically appropriate for their analysis. The objective of the present study was to compare models for identification of QTL associated with survival traits, in particular when combined with selective genotyping. Data were simulated to model the survival distribution of a population of chickens challenged with Marek disease virus. Cox proportional hazards (CPH), linear regression (LR), and Weibull models were compared for their appropriateness to analyze the data, ability to identify associations of marker alleles with survival, and estimation of effects when all individuals were genotyped (full genotyping) and when selective genotyping was used. Little difference in power was found between the CPH and the LR model for low censoring cases for both full and selective genotyping. The simulated data were not transformed to follow a Weibull distribution and, as a result, the Weibull model generally resulted in less power than the other two models and overestimated effects. Effect estimates from LR and CPH were unbiased when all individuals were genotyped, but overestimated when selective genotyping was used. Thus, LR is preferred for analyzing survival data when the amount of censoring is low because of ease of implementation and interpretation. Including phenotypic data of non-genotyped individuals in selective genotyping analysis increased power, but resulted in LR having an inflated false positive rate, and therefore the CPH model is preferred for this scenario, although transformation of the data may also make the Weibull model appropriate for this case. The results from the research presented herein are directly applicable to interval mapping analyses.     

Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6

Background

In pig, a number of experiments have been set up to identify QTL and a multitude of chromosomal regions harbouring genes influencing traits of interest have been identified. However, the mapping resolution remains limited in most cases and the detected QTL are rather inaccurately located. Mapping accuracy can be improved by increasing the number of phenotyped and genotyped individuals and/or the number of informative markers. An alternative approach to overcome the limited power of individual studies is to combine data from two or more independent designs.

Methods

In the present study we report a combined analysis of two independent design (a French and a Dutch F2 experimental designs), with 2000 F2 individuals. The purpose was to further map QTL for growth and fatness on pig chromosomes 2, 4 and 6. Using QTL-map software, uni- and multiple-QTL detection analyses were applied separately on the two pedigrees and then on the combination of the two pedigrees.

Results

Joint analyses of the combined pedigree provided (1) greater significance of shared QTL, (2) exclusion of false suggestive QTL and (3) greater mapping precision for shared QTL.

Conclusions

Combining two Meishan x European breeds F2 pedigrees improved the mapping of QTL compared to analysing pedigrees separately. Our work was facilitated by the access to raw phenotypic data and DNA of animals from both pedigrees and the combination of the two designs with the addition of new markers allowed us to fine map QTL without phenotyping additional animals.     

A gene frequency model for QTL mapping using Bayesian inference

Background

Information for mapping of quantitative trait loci (QTL) comes from two sources: linkage disequilibrium (non-random association of allele states) and cosegregation (non-random association of allele origin). Information from LD can be captured by modeling conditional means and variances at the QTL given marker information. Similarly, information from cosegregation can be captured by modeling conditional covariances. Here, we consider a Bayesian model based on gene frequency (BGF) where both conditional means and variances are modeled as a function of the conditional gene frequencies at the QTL. The parameters in this model include these gene frequencies, additive effect of the QTL, its location, and the residual variance. Bayesian methodology was used to estimate these parameters. The priors used were: logit-normal for gene frequencies, normal for the additive effect, uniform for location, and inverse chi-square for the residual variance. Computer simulation was used to compare the power to detect and accuracy to map QTL by this method with those from least squares analysis using a regression model (LSR).

Results

To simplify the analysis, data from unrelated individuals in a purebred population were simulated, where only LD information contributes to map the QTL. LD was simulated in a chromosomal segment of 1 cM with one QTL by random mating in a population of size 500 for 1000 generations and in a population of size 100 for 50 generations. The comparison was studied under a range of conditions, which included SNP density of 0.1, 0.05 or 0.02 cM, sample size of 500 or 1000, and phenotypic variance explained by QTL of 2 or 5%. Both 1 and 2-SNP models were considered. Power to detect the QTL for the BGF, ranged from 0.4 to 0.99, and close or equal to the power of the regression using least squares (LSR). Precision to map QTL position of BGF, quantified by the mean absolute error, ranged from 0.11 to 0.21 cM for BGF, and was better than the precision of LSR, which ranged from 0.12 to 0.25 cM.

Conclusions

In conclusion given a high SNP density, the gene frequency model can be used to map QTL with considerable accuracy even within a 1 cM region.     

Gene expression profiling of early intervertebral disc degeneration reveals a down-regulation of canonical Wnt signaling and caveolin-1 expression: implications for development of regenerative strategies

Introduction

Early degeneration of the intervertebral disc (IVD) involves a change in cellular differentiation from notochordal cells (NCs) in the nucleus pulposus (NP) to chondrocyte-like cells (CLCs). The purpose of this study was to investigate the gene expression profiles involved in this process using NP tissue from non-chondrodystrophic and chondrodystrophic dogs, a species with naturally occurring IVD degeneration.

Methods

Dual channel DNA microarrays were used to compare 1) healthy NP tissue containing only NCs (NC-rich), 2) NP tissue with a mixed population of NCs and CLCs (Mixed), and 3) NP tissue containing solely CLCs (CLC-rich) in both non-chondrodystrophic and chondrodystrophic dogs. Based on previous reports and the findings of the microarray analyses, canonical Wnt signaling was further evaluated using qPCR of relevant Wnt target genes. We hypothesized that caveolin-1, a regulator of Wnt signaling that showed significant changes in gene expression in the microarray analyses, played a significant role in early IVD degeneration. Caveolin-1 expression was investigated in IVD tissue sections and in cultured NCs. To investigate the significance of Caveolin-1 in IVD health and degeneration, the NP of 3-month-old Caveolin-1 knock-out mice was histopathologically evaluated and compared with the NP of wild-type mice of the same age.

Results

Early IVD degeneration involved significant changes in numerous pathways, including Wnt/β-catenin signaling. With regard to Wnt/β-catenin signaling, axin2 gene expression was significantly higher in chondrodystrophic dogs compared with non-chondrodystrophic dogs. IVD degeneration involved significant down-regulation of axin2 gene expression. IVD degeneration involved significant down-regulation in Caveolin-1 gene and protein expression. NCs showed abundant caveolin-1 expression in vivo and in vitro, whereas CLCs did not. The NP of wild-type mice was rich in viable NCs, whereas the NP of Caveolin-1 knock-out mice contained chondroid-like matrix with mainly apoptotic, small, rounded cells.

Conclusions

Early IVD degeneration involves down-regulation of canonical Wnt signaling and Caveolin-1 expression, which appears to be essential to the physiology and preservation of NCs. Therefore, Caveolin-1 may be regarded an exciting target for developing strategies for IVD regeneration.     

SIMPLE: implementation of recommendations from international evidence-based guidelines on caesarean sections in the Netherlands. Protocol for a controlled before and after study

Background

Caesarean section (CS) rates are rising worldwide. In the Netherlands, the most significant rise is observed in healthy women with a singleton in vertex position between 37 and 42 weeks gestation, whereas it is doubtful whether an improved outcome for the mother or her child was obtained. It can be hypothesized that evidence-based guidelines on CS are not implemented sufficiently. Therefore, the present study has the following objectives: to develop quality indicators on the decision to perform a CS based on key recommendations from national and international guidelines; to use the quality indicators in order to gain insight into actual adherence of Dutch gynaecologists to guideline recommendations on the performance of a CS; to explore barriers and facilitators that have a direct effect on guideline application regarding CS; and to develop, execute, and evaluate a strategy in order to reduce the CS incidence for a similar neonatal outcome (based on the information gathered in the second and third objectives).

Methods

An independent expert panel of Dutch gynaecologists and midwives will develop a set of quality indicators on the decision to perform a CS. These indicators will be used to measure current care in 20 hospitals with a population of 1,000 women who delivered by CS, and a random selection of 1,000 women who delivered vaginally in the same period. Furthermore, by interviewing healthcare professionals and patients, the barriers and facilitators that may influence the decision to perform a CS will be measured. Based on the results, a tailor-made implementation strategy will be developed and tested in a controlled before-and-after study in 12 hospitals (six intervention, six control hospitals) with regard to effectiveness, experiences, and costs.

Discussion

This study will offer insight into the current CS care and into the hindering and facilitating factors influencing obstetrical policy on CS. Furthermore, it will allow definition of patient categories or situations in which a tailor-made implementation strategy will most likely be meaningful and cost effective, without negatively affecting the outcome for mother and child.

Trial registration

http://www.clinicaltrials.gov: NCT01261676     

Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees

Background

In the pig, multiple QTL associated with growth and fatness traits have been mapped to chromosome 2 (SSC2) and among these, at least one shows paternal expression due to the IGF2-intron3-G3072A substitution. Previously published results on the position and imprinting status of this QTL disagree between analyses from French and Dutch F2 crossbred pig populations obtained with the same breeds (Meishan crossed with Large White or Landrace).

Methods

To study the role of paternal and maternal alleles at the IGF2 locus and to test the hypothesis of a second QTL affecting backfat thickness on the short arm of SSC2 (SSC2p), a QTL mapping analysis was carried out on a combined pedigree including both the French and Dutch F2 populations, on the progeny of F1 males that were heterozygous (A/G) and homozygous (G/G) at the IGF2 locus. Simulations were performed to clarify the relations between the two QTL and to understand to what extent they can explain the discrepancies previously reported.

Results

The QTL analyses showed the segregation of at least two QTL on chromosome 2 in both pedigrees, i.e. the IGF2 locus and a second QTL segregating at least in the G/G F1 males and located between positions 30 and 51 cM. Statistical analyses highlighted that the maternally inherited allele at the IGF2 locus had a significant effect but simulation studies showed that this is probably a spurious effect due to the segregation of the second QTL.

Conclusions

Our results show that two QTL on SSC2p affect backfat thickness. Differences in the pedigree structures and in the number of heterozygous females at the IGF2 locus result in different imprinting statuses in the two pedigrees studied. The spurious effect observed when a maternally allele is present at the IGF2 locus, is in fact due to the presence of a second closely located QTL. This work confirms that pig chromosome 2 is a major region associated with fattening traits.     

A polycystin‐2 (TRPP2) dimerization domain essential for the function of heteromeric polycystin complexes

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in two genes, PKD1 and PKD2, which encode polycystin‐1 (PC1) and polycystin‐2 (PC2), respectively. Earlier work has shown that PC1 and PC2 assemble into a polycystin complex implicated in kidney morphogenesis. PC2 also assembles into homomers of uncertain functional significance. However, little is known about the molecular mechanisms that direct polycystin complex assembly and specify its functions. We have identified a coiled coil in the C‐terminus of PC2 that functions as a homodimerization domain essential for PC1 binding but not for its self‐oligomerization. Dimerization‐defective PC2 mutants were unable to reconstitute PC1/PC2 complexes either at the plasma membrane (PM) or at PM‐endoplasmic reticulum (ER) junctions but could still function as ER Ca²⁺‐release channels. Expression of dimerization‐defective PC2 mutants in zebrafish resulted in a cystic phenotype but had lesser effects on organ laterality. We conclude that C‐terminal dimerization of PC2 specifies the formation of polycystin complexes but not formation of ER‐localized PC2 channels. Mutations that affect PC2 C‐terminal homo‐ and heteromerization are the likely molecular basis of cyst formation in ADPKD.     

Genome-wide association study of insect bite hypersensitivity in two horse populations in the Netherlands

Background

Insect bite hypersensitivity is a common allergic disease in horse populations worldwide. Insect bite hypersensitivity is affected by both environmental and genetic factors. However, little is known about genes contributing to the genetic variance associated with insect bite hypersensitivity. Therefore, the aim of our study was to identify and quantify genomic associations with insect bite hypersensitivity in Shetland pony mares and Icelandic horses in the Netherlands.

Methods

Data on 200 Shetland pony mares and 146 Icelandic horses were collected according to a matched case–control design. Cases and controls were matched on various factors (e.g. region, sire) to minimize effects of population stratification. Breed-specific genome-wide association studies were performed using 70 k single nucleotide polymorphisms genotypes. Bayesian variable selection method Bayes-C with a threshold model implemented in GenSel software was applied. A 1 Mb non-overlapping window approach that accumulated contributions of adjacent single nucleotide polymorphisms was used to identify associated genomic regions.

Results

The percentage of variance explained by all single nucleotide polymorphisms was 13% in Shetland pony mares and 28% in Icelandic horses. The 20 non-overlapping windows explaining the largest percentages of genetic variance were found on nine chromosomes in Shetland pony mares and on 14 chromosomes in Icelandic horses. Overlap in identified associated genomic regions between breeds would suggest interesting candidate regions to follow-up on. Such regions common to both breeds (within 15 Mb) were found on chromosomes 3, 7, 11, 20 and 23. Positional candidate genes within 2 Mb from the associated windows were identified on chromosome 20 in both breeds. Candidate genes are within the equine lymphocyte antigen class II region, which evokes an immune response by recognizing many foreign molecules.

Conclusions

The genome-wide association study identified several genomic regions associated with insect bite hypersensitivity in Shetland pony mares and Icelandic horses. On chromosome 20, associated genomic regions in both breeds were within 2 Mb from the equine lymphocyte antigen class II region. Increased knowledge on insect bite hypersensitivity associated genes will contribute to our understanding of its biology, enabling more efficient selection, therapy and prevention to decrease insect bite hypersensitivity prevalence.     

Low sclerostin levels: a predictive marker of persistent inflammation in ankylosing spondylitis during anti-tumor necrosis factor therapy?

Introduction

Sclerostin levels have been reported to be low in ankylosing spondylitis (AS), but there is no data regarding the possible role of this Wnt inhibitor during anti-tumor necrosis factor (TNF) therapy. The present study longitudinally evaluated sclerostin levels, inflammatory markers and bone mineral density (BMD) in AS patients under anti-TNF therapy.

Methods

Thirty active AS patients were assessed at baseline, 6 and 12 months after anti-TNF therapy regarding clinical parameters, inflammatory markers, BMD and baseline radiographic damage (mSASSS). Thirty age- and sex-matched healthy individuals comprised the control group. Patients'' sclerostin levels, sclerostin binding low-density lipoprotein receptor-related protein 6 (LRP6) and BMD were evaluated at the same time points and compared to controls.

Results

At baseline, AS patients had lower sclerostin levels (60.5 ± 32.7 vs. 96.7 ± 52.9 pmol/L, P = 0.002) and comparable sclerostin binding to LRP6 (P = 0.387) than controls. Improvement of Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), Bath Ankylosing Spondylitis Functional Index (BASFI), Bath Ankylosing Spondylitis Metrology Index (BASMI), Ankylosing Spondylitis quality of life (ASQoL) was observed at baseline vs. 6 vs. 12 months (P < 0.01). Concomitantly, a gradual increase in spine BMD (P < 0.001) and a positive correlation between baseline mSASSS and spine BMD was found (r = 0.468, P < 0.01). Inflammatory parameters reduction was observed comparing baseline vs. 6 vs. 12 months (P <0.01). Sclerostin levels progressively increased [baseline (60.5 ± 32.7) vs. 6 months (67.1 ± 31.9) vs. 12 months (72.7 ± 32.3) pmol/L, P <0.001]. At 12 months, the sclerostin levels remained significantly lower in patients compared to controls (72.7 ± 32.3 vs. 96.70 ± 52.85 pmol/L, P = 0.038). Moreover, sclerostin serum levels at 12 months were lower in the 10 patients with high C reactive protein (CRP) (≥ 5 mg/l) compared to the other 20 patients with normal CRP (P = 0.004). Of note, these 10 patients with persistent inflammation also had lower sclerostin serum levels at baseline compared to the other patients (P = 0.023). Univariate logistic regression analysis demonstrated that AS patients with lower sclerostin serum levels had an increased risk to have high CRP at 12 months (odds ratio = 7.43, 95% CI 1.23 to 45.01, P = 0.020) than those with higher sclerostin values.

Conclusions

Persistent low sclerostin levels may underlie continuous inflammation in AS patients under anti-TNF therapy.     

Tattoo pigment in an axillary lymph node simulating metastatic malignant melanoma

We report a case of axillary lymphadenopathy thirty years after a decorative tattoo clinically mimicking metastatic melanoma. The importance of relying on histological confirmation of metastatic disease before altering extent of surgery is discussed. The importance of recording presence of decorative tattoos is stressed.